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errors to produce command signals representing a force 
vector F and a torque vector T, each having three or- 
thogonal components, which represent the forces and 
torques which are to be generated by the thrusters. The 
thrusters may include magnetic torquer or reaction 
wheels. Six difference equations are generated, three 
having the form 

N 

A,* = F\ - Sj bip.fi 
and 

A 7 

A/4-3 — Ti — OifLj- 

where ay is the maximum torque which the } th thruster 
can produce, by is the maximum force which the } ih 
thruster can produce, and ay is a variable representing 
the throttling factor of the } th thruster, which may range 
from zero to unity. The six equations are summed to 
produce a single scalar equation relating variables ay to 
a performance index Z: 
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[57] ABSTRACT 

A spacecraft attitude and/or velocity control system 
includes a controller which responds to at least attitude 


6 N 

Z = - 2 | A/| - 2 a/ 

i=l /= 1 

Those values of a which maximize the value of Z are 
determined by a method for solving linear equations, 
such as a linear programming method. The Simplex 
method may be used. The values of ay are applied to 
control the corresponding thrusters. 
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SPACECRAFT ATTITUDE AND VELOCITY 
CONTROL SYSTEM 

The invention described herein was made in the per- 5 
formance of work under NASA Contract No. NAS5- 
32000 and is subject to the provisions of Section 305 of 
the National Aeronautics and Space Act of 1958, as 
amended (72 Stat 435; 42 USC 2457). 

BACKGROUND OF THE INVENTION 

This invention relates to the control of the thrusters 
of a spacecraft, and more particularly, to attitude and 
velocity control by the use of throttleable or pulse- 
modulated thrusters. 

Spacecraft are widely used for communications, 
earth sensing and exploration, vehicle locating, and for 
surveillance. All of these uses require that the orienta- 
tion of the spacecraft in space, and possibly its location, jq 
be accurately controlled. The spacecraft attitude may 
be controlled by magnetic torquers, by momentum or 
reaction wheels, or by the use of thrusters. 

Since each spacecraft at launch has a slightly differ- 
ent configuration than other spacecraft, its thrusters 25 
will be placed at slightly different locations relative to 
the center of gravity than in other spacecraft. Also, the 
specific impulse characteristics of the thrusters may 
differ from one to the next. As a consequence, the logic 
for controlling the thruster system is ordinarily custom- 30 
written for each spacecraft. During the construction 
phase for a spacecraft, it may be found to be necessary 
to alter the configuration or thruster characteristics, 
which may require costly rewriting of software. 

Control arrangements ordinarily provide torques by 35 
preselected combinations of thrusters. The selected 
thrusters are ordinarily located on opposite sides of the 
center of gravity. When in orbit, a spacecraft thruster 
may fail or may have a significant change in specific 
impulse characteristics. When a thruster fails, those 40 
combinations of torque for which the failed thruster is 
an element cannot be used. This may adversely affect or 
limit control of the spacecraft. 

A spacecraft velocity and attitude control arrange- 45 
ment is desired in which a change in specific impulse of 
a particular thruster need not necessarily eliminate its 
use in controlling the spacecraft. 

SUMMARY OF THE INVENTION 

50 

A method for controlling the attitude and velocity of 
a spacecraft includes the steps of generating sensed 
signals representative of pitch, roll and yaw. A velocity 
error signal may also be generated. Command signals 
are generated representing force and torque vectors F 55 
and T, respectively, which are to be generated by the 
thrusters. The force vector includes mutually orthogo- 
nal components Fj, F 2 and F 3 , and the torque vector 
includes mutually orthogonal components Ti, T 2 and 
T 3 . When unidimensional corrections are to be made, 60 
some or most of these components may take on zero 
magnitude. Six difference equations A of the form: 

N 

A/ ~ F\ - ^ X ^ bijaf, 65 


2 

-continued 

N 

A+3 = Ti - Qifij, 

are generated, where: a and b are the maximum torque 
and force, respectively, which the j r/> thruster can pro- 
duce, and a is the throttle setting of the j th thruster, 
which may take on values ranging from zero to unity. 
The difference signals A are summed to form a single 
scalar equation relating variable a to a scalar perfor- 
mance index Z. The single scalar equation is solved for 
that value of ay which maximizes Z, and the thrusters are 
controlled in a manner directly related to the corre- 
sponding value of ay. 

DESCRIPTION OF THE DRAWING 

FIG. la is a perspective or isometric view of a space- 
craft including a plurality of thrusters, FIG. lb illus- 
trates a simplified side view of the spacecraft of FIG. 
la, and FIG. 1c is a simplified side view of the space- 
craft of FIG. la; and 

FIG. 2 is a flow chart illustrating a control method 
according to the invention. 

DESCRIPTION OF THE INVENTION 

FIG. la illustrates a spacecraft which is controlled in 
accordance with the invention. In FIG. la, a spacecraft 
10 having a body 12 includes a omnidirectional antenna 
illustrated as 14 which may be used for communication 
and command, and also includes attitude sensors for 
roll, pitch and yaw, which are illustrated together as a 
block 16. Roll, pitch and yaw axes 4, 6 and 8 are associ- 
ated with spacecraft body 12. Attitude signals represen- 
tative of roll, pitch and yaw are coupled (by means 
which are not illustrated) to a controller illustrated as a 
block 18 for generating thruster control signals. A plu- 
rality of thrusters illustrated as 20a-20r are associated 
with spacecraft body 12 and are connected (by means 
not illustrated) to controller 18 for being controlled 
thereby. 

FIG. lb illustrates spacecraft 10 in simplified form, 
looking along the yaw axis. Elements of FIG. lb corre- 
sponding to those of FIG. la are designated by the same 
reference numerals. FIG. 16 is simplified, in that only 
thrusters 20/ and 20/ of FIG. la are illustrated, together 
with additional thrusters 20s, 20/ which are not visible 
in FIG. la. Thrusters 20/, 20), 20s and 20/ when oper- 
ated together are capable of generating a torque in the 
direction of arrow 22 about a center of gravity (CG) 24. 

FIG. lc is similar to FIG. lb, and illustrates additional 
thrusters 20 u, 20 v which, together with thrusters 20s, 
20/ are capable of creating a force for accelerating the 
spacecraft in the direction of arrow 26. Since the center 
of gravity 24 is not equidistant between the projections 
of the direction of thrust of thrusters 20s, 20/, 20 u and 
20 v, a torque will be created. Such a torque may be 
undesirable when only acceleration in a given direction 
is desired. Such a torque may be countered by operating 
thrusters 20/i and 20 k sufficiently to oppose the torque. 
Alternatively, translation may be accomplished without 
the torque by producing less thrust from thrusters 20s, 
20/ than from 20a, 20v. 

The attitude of the spacecraft can be determined in 
any conventional manner. A method for determining 
the attitude of a spacecraft using a polar star sensor is 
described in U.S. Pat. application Ser. No. 07/511,169, 
filed Apr. 19, 1990 in the name of Paluszek. A method 


and 
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for determining the attitude of a spacecraft by the use of 
Doppler information measured from a ground station is 
described in U.S. Pat application Ser. No. 07/397,939, 
filed Aug. 24, 1989 now U.S. Pat. No. 5,070,338 in the 
name of Cohen. Spacecraft velocity may also be deter- 5 
mined in any conventional manner, such as by the use of 
radar either from the spacecraft or from a ground sta- 
tion. If a ground station determines the velocity or the 
attitude, the information may be transmitted to the 
spacecraft for reception by antenna 12 and for commu- io 
nication to controller 18. In general, attitude control 
requires that the spacecraft be torqued, and station 
keeping requires that forces be imported to the space- 
craft. As mentioned in relation to FIGS, lb and lc, 
application of forces in turn may require application of 15 
torque to correct for residual torques imparted by the 
thrusters. 

Controller 18 of FIG. la responds to the attitude 
signals and, as appropriate, to the velocity signals in a 
conventional manner to produce command signals rep- 20 
resenting a force vector F and a torque vector T, which 
represent the force and torque vectors which are to be 
generated by the thrusters. The force vector F includes 
mutually orthogonal components Fi, F2 and F3. The 
torque vector T includes mutually orthogonal torque 25 
components Ti, T2and T3. Naturally, depending upon 
the attitude error, some or most of the components may 
take on zero value for all or part of a maneuver. The 
method illustrated in FIG. 2 shows the reading of the 
attitude control signals as a block 210, and the genera- 30 
tion of vectors F and T as a block 212. 

In accordance with the invention, the control method 
includes the generation of six different equations A: 


N 

(I) 35 

A\ - F\ - j2 j b\ja.j 

N 

A 2 = Fi - ,2^ b2jaj 

(2) 


40 

N 

(3) 

A3 = ^3 - 2j 63/0/ 

N 

A 4 - 21 - l/tty 

(4) 

45 

N 

As = Ti - 2^ a2jaj 

(5) 

N 

(6) 

A6 - 73 - 2 ^ azjCLj 
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single scalar equation relating variable a to a scalar 
performance index Z. 


Block 220 of FIG. 2 represents the solving of equa- 
tion 7 for those values of a which maximize Z. The 
solution is achieved by some method for solving linear 
equations, such as a linear programming method. In 
particular, the Simplex method may be used. The Sim- 
plex method is described in standard texts, such as “Nu- 
merical Recipes,” William H. Press et al., Cambridge 
University Press, 1986, pp. 289-293. 

Finally, as illustrated by block 222 of FIG. 2, the 
values of a as determined in block 216 are applied for 
control of the corresponding thrusters to achieve the 
desired control result. This requires only a proportional 
throttling of the thrusters. If pulse-width modulation 
control is used, the ON to OFF durations are set to the 
required proportion for each thruster. 

For North-South station keeping maneuvers for a 
geosynchronous satellite, to correct orbital inclination 
errors, translational forces in a North-South direction 
may be necessary. If the spacecraft is too far North, 
maximum AV in a Southerly direction may be neces- 
sary. In principle, this requires that the thrusters on the 
North face of the spacecraft be operated to provide 
thrust along the pitch axis. However, differences in 
thrust among the North-face thrusters, attitude errors 
and the like may result in net thrust in undesired direc- 
tions. Prior to starting the maneuver, attitude-sensing 
gyros are initialized by the use of Earth sensor data. 
Attitude errors sensed by the gyros during the maneu- 
ver are applied to the control arrangement, correspond- 
ing to block 212 of FIG. 2. The control arrangement 
may use a proportional-derivative algorithm, or a pro- 
portional-integral-derivative (PID) algorithm. A PID 
controller is described in U.S. Pat. application Ser. No. 
07/459,627, filed Jan. 2, 1990 now U.S. Pat. No. 
5,025,381 in the name of Goodzeit et al. The controller 
calculates control torques T about pitch, roll and yaw 
axes which are required for stabilizing the spacecraft. 
Torques T are used in block 214 to generate A4, A5 and 
A6- For this particular application, the cost functional Z 
produced in block 218 of FIG. 1 has the form: 

6 N 2 N (*) 

Z - - 2 | A/| -f 2 aj- 2 aj 
*= 4 J^P 1 


In equations 1-6, ay is the maximum torque which the 
thruster can produce, and relates to the force which 
the j lk thruster can produce and its distance from the 
center of gravity. There are a total of N thrusters. The 55 
term by is the maximum force which the j th thruster can 
produce. The values of ay and by are predetermined, and 
may be established by prelaunch measurements or by 
on-orbit tests performed from time to time. Variable ay 
represents the throttle setting of the j th thruster, which 60 
may take on values ranging from zero to unity. The step 
of forming equations 1-6 is illustrated in FIG. 2 as a 
block 214. The step represented by block 214 uses val- 
ues a / and b / which are read from memory, as repre- 
sented by block 216. 65 

Block 218 of FIG. 2 illustrates a further step of the 
control method according to the invention, namely that 
of summing the difference equations 1-6 to produce a 


which is generated by substituting 
*2 3 

2 ay for 2 |A/| 

J~P < = 1 

in equation 6. 

In equation 7, thrusters with indices ranging from p to 
N2 are oriented for thrust in the direction opposite to 
the desired AV direction. The values of a are applied to 
the corresponding thrusters, so that attitude control 
about roll and yaw axes is accomplished by throttling or 
modulating the North-face thrusters. Attitude control 
about the pitch axis (the AV axis) is exerted by pulsing 
ON the East and West face thrusters. 

Thus, maximizing the cost functional Z will maximize 
AV, and will minimize errors between the desired con- 
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trol torques and the actual control torques. The cost 
functional equation (equation 6 or 7) is solved during 
each control cycle using linear programming method to 
maximize Z, and the throttle settings ay always remain 
in the range of 0 to 1 . The throttle values ay are used to 5 
set the pulse width (or other throttle control) for each 
thruster, so that the pulse width equals ray, where t is 
the control period. 

The method according to the invention has a further io 
advantage, in that the control algorithm which gener- 
ates the desired F and T can be common among a plu- 
rality of spacecraft, thereby achieving a cost reduction 
and an increase in effective reliability. Since only the 
thruster characteristics and their number and location 15 
relative to the center of gravity change from spacecraft 
to spacecraft, only the values of a y and by need be 
changed in the equations represented by block 216 of 
FIG. 2, and the control logic may be reused but for 
these minor changes. 

Other embodiments of the invention will be apparent 
to those skilled in the art. In particular, this approach 
could be used with other actuators such as magnetic 
torquers and reaction wheels. It could be used with 25 
combinations of momentum exchange devices and low 
thrust thrusters. It could be used solely as part of a 
momentum unloading system in which momentum ex- 
change devices were used for control. In addition, pre- 
defined weighting factors could be added to increase or 
decrease the significance of a term in the cost equation. 
What is claimed is: 

1. A method for controlling the attitude and velocity 
of a spacecraft, comprising the steps of: 
generating sensed signals representative of pitch, roll, 
yaw and at least velocity error of said spacecraft; 
from said sensed signals, generating command signals 
representing a force vector F, and torque vector T 
to be generated by said thrusters, said force vector 40 
F including mutually orthogonal components Fj, 

F 2 and F 3 , and said torque vector T including mu- 
tually orthogonal components Ti, T 2 and T 3 , some 
or all of which said Fi, F 2 , Fjand T\, T 3 and T 3 can 
from time to time take on zero magnitude; 
generating difference signals representing six differ- 
ence equations A of the form: 


30 


-continued 

N 

A6 - 73 - 2 aijcLj 


where: 

N is the total number of thrusters; 

aiy, a 2 yand a 3 yare the maximum torques which the j th 
thrusters can produce about the 1st, 2nd and 3nd 
axes; 

byb 2 yb 3 yare the maximum forces which the j th thrust- 
ers can produce along the 1 st, 2nd and 3rd axes; 
and 

aj is a variable representing the throttle setting of the 
j th thruster, which can take on volumes ranging 
from zero to unity; 

summing said difference signals to form a single sca- 
lar equation relating the variable aj to a scalar 
performance index Z; 

solving said scalar equation for a value of the variable 
aj which maximizes said Z; and 

controlling each thruster in a manner directly related 
to the corresponding value of said aj. 

2. A method according to claim 1 wherein said solv- 
ing step comprises the step of applying a linear pro- 
gramming method algorithm for solving said scalar 
equation. 

3. A method according to claim 2 wherein said linear 
programming method algorithm comprises a Simplex 
algorithm. 

4. A method according to claim 1, wherein said sum- 
ming step forms a single scalar equation in which the 
equation defining said scalar Z includes at least the term 


35 


6 

-1 | A/| 
i=m 


N 2 

+ 2 ay 
j-P 


45 


where thrusters which indices ranging from P to N 2 
are oriented for thrust in the direction opposite to 
the desired AV direction, and m=l. 

5. A method according to claim 4, wherein p = 1 and 
N 2 equals the total number of thrusters. 

6 . A method according to claim 4, wherein said equa- 
tion defining said scalar Z includes, in addition to the 
term 


N 

Ai = Fi — b\jcLj 

N 

a 2 = f 2 - Hp-j 

N 

A3 = F 2 - 2^ bsjaj 
N 

A4 = Ti - fll/xy 

N 

A 5 = r 2 -^2^ a 2 jCLj 
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6 

-.2 | A/| 

t—m 


N 2 

+ 2 < 
J=P 


the additional term 


55 


60 


N\ 

— 2 ay 
7=1 1 

where p=Ni, and thrusters having indices ranging 
from 1 to Nj do not have their thrust vectors 
aligned with a desired direction of acceleration, 
while thrusters having indices ranging from Ni to 
N 2 have their thrust vectors aligned with the de- 
sired direction of acceleration. 
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